Análise Matemática IV 

Problemas para as Aulas Práticas 



Semana 8 



1. Para cada uma das seguintes equações diferenciais, esboce o campo de direcções e trace 
os respectivos tipos de soluções . 



(a) y' = 



ty 



i + t 2 ' 

(c)y' = y(l-y 2 ) 



(b)y'=(2-y)(y-l), 



( d ) y = 



y-t 1 



(e) y' = sen(y - t) 



Resolução: 

ty 



(a) y' = 



i + t 2 




(b)y'=(2-y)(y-l) 



(c)y' = y(l-y 2 ) 



(d) y' = 



y + t 
y-t 
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(e) y = sen(y - t) 




2. Mostre que existe uma solução de classe C 1 para o problema de valor inicial 



' dy 



= QtUy 



?/„,2 



dt 
,2/(0) = 0, 

diferente da solução y(t) = 0, Vi £ K. Explique porque é que isto não contradiz o teorema de 
Picard. 

Resolução: 

Como referido no enunciado, uma das soluções do PVI é a solução constante y{t) = 0. Por 
outro lado, se y(t) ^ 0, a equação pode ser escrita na forma 



-2/3 d V n, „ d 

y ' — - = bt <í=> — 

dt dt 



( í ' y- 2/3 dy) = 6í & 3y 1/s = 3t 2 + c & y(t) = {t 2 + kf 



Esta solução polinomial, logo é diferenciável em KL A condição inicial y(0) = é satisfeita se 
tomarmos k = 0, logo outra solução do PVI é 

y{t) = t\ 

Para verificar que não há contradição com o Teorema de Picard, note-se que, sendo f(t,y) = 
ôtyy 2 , f é contínua em M 2 e 



lim — — = lim Aty' 



-1/3 _ 



00 



y^o dy y^o 

pelo que é de esperar que f(t,y) não seja localmente lipschitziana em ordem a y em qualquer 
subconjunto compacto de M 2 que contenha (ío,2/o) = (0,0). De facto, se \t\ < a e \y\ < 0, 



\f(t,y)- f(t,x)\ = \6t^-6t\/^\ = |6t||y 2/3 -x 2/3 | = 6|í| 
e é fácil de verificar que para y, x numa vizinhança de a função 

y 2/3 _ x 2/3 



V 



2/3 _ ^2/3 



y- x 



\y-x\ 



y-x 



não é limitada. Concluímos então que a continuidade de / implica existência de solução do 
PVI, mas o facto de não ser localmente lipschitziana em ordem a y numa vizinhança de (0, 0) 
não assegura a unicidade de solução do PVI. 



3. Mostre que o problema de valor inicial 

dy . 

Tt = y/ 

,y(o) = o, 

tem infinitas soluções, e explique porque esse facto não contradiz o Teorema de Picard. 

Resolução: 

Começamos por verificar que a solução constante y(t) = é solução do PVI. Por outro lado, 
se y(t) ^0a equação pode ser escrita na forma 



V 



-1/2^2/ = 1 ^ £ 

dt dt 



(Jy- 1/2 dy} =1 & 2y 1 / 2 = t + c O y(t)= (t + k 



Visto termos obtido uma função polinomial verifica-se que y é diferenciável em M, e y(0) = 
implica que outra solução do PVI é 

yit) = - 

Podemos agora utilizar "cortar" e "colar" entre estas duas soluções para criar novas soluções 
do PVI. Isto é, para t > 0, defina-se 







se t < t 



o 



yt (t) = 



,t ío\ 2 

se t > tn 

v 2 2/ 



Verifica-se que y to é diferenciável em M, verifica a equação diferencial (dado que e 
verificam) e yt o (0) = 0. De igual modo, para cada sq < 

se t > sq 

y so (t) = 



t t \2 



{---) set< s 



é também solução do PVI. 

Finalmente, o facto de existirem uma infinidade de soluções deve-se a que a função f(t, y) = 
^/y é contínua em y > 0, mas não é localmente lipschitziana em ordem a y em qualquer conjunto 
compacto de M? que contenha a origem (0, 0). De facto, temos que: 



\f(t,x)-f(t,y)\ = 



\fx~ y/V 



x-y 



\x-y\, 



onde a função 



Vx- y/y 



x-y 



não +e limitada para x, y num vizinhança qualquer da origem. 



4. Considere o seguinte problema de valor inicial 

'1 +\ dy 1 2 
l-t)y- = l-y 

y(l/2) = 2, 

(i) Determine uma solução do PVI, e justifique que essa é a única solução do problema 
definida para t numa vizinhança de 1/2. 

(ii) Mostre que o PVI admite um número infinito de soluções definidas em M. 

(iii) Diga, justificando, porque não há contradição ao Teorema de Picard. 

Resolução: 

(i) Começamos por observar que a equação diferencial faz sentido para qualquer í Ele 
qualquer t/el. Trata-se de uma equação separável, pelo que para í / 1 ey 2 ^ 1: 



y 2 -ldí í-1 dtj \y 2 -l *J í-1 

O -log(y 2 -l) = log(í-l) + c 

O y 2 = fc(í - l) 2 + 1 

O y(í) = v 7 ^" 1) 2 + 1 ou y(í) = -^k(t- l) 2 + 1 
Dado que y(l/2) = 2 > 0, a solução do PVI é 

y(t) = ^l + 12(t-iy (1) 

Para mostrar que é a única solução do PVI teremos que verificar que 

i-y 2 



f(t,y) = 



y(i-*) 



verifica as condições do Teorema de Picard em certo conjunto D <Z M. 2 contendo a condição 
inicial (to,yo) = (1/2,2). O domínio de / é D = {(t,y) : ?/ / e í ^ 1} e é óbvio que 
(to,Vo) = (1/2) 2) G D. Por outro lado dado que / é uma função racional, em D, f é contínua 

d_£ 1 + y 2 

dy [hy) y 2 (l-t) 

é também contínua em D (pelo que / é localmente lipschitziana em relação a y em D). Estamos 
então nas condições do Teorema de Picard e concluir que o PVI admite solução única numa 
vizinhança de t = 1/2. Pela unicidade conclui-se que a solução tem que ser dada por (1). 

(ii) Comecemos por calcular o intervalo máximo da solução calculada na alínea anterior. 
Sabemos que I =]a, /3[ em que 

to = 1/2 G /; 

quando t —>■ a + ou (t,y(t)) —> dD ou \y(t)\ —> oo; 



quando t —> (3 ou (t,y(t)) —> dD ou \y(t)\ —> oo. 

Visto o domínio de diferenciabilidade de y(t) ser IR, podemos desde já concluir que y(t) não 
explode em tempo finito. Por isso, o único problema que pode surgir é o de (t,y(t)) atingir a 
fronteira de D, isto é, quando t = 1 ou y(t) = 0. Mais uma vez, pela expressão de y(t) podemos 
concluir que y(t) ^ para todo íel, pelo que o único problema é mesmo t = 1. Tem-se então 
que o intervalo máximo de solução única é 

I=]-oo,l[ 

e 

y(t) —> 1 quando t —>■ 1~ 

Por outro lado, sabemos que 

y k (t) = y/k(t -1) 2 + 1 

é solução da equação diferencial e 

y k il) = 1 , VA: > 

Assim, yk{t) é solução do problema 



(l-t)y^- = l-y 

v >y dt 



2 



lim i _ >1 + y(t) = 1 ou i/(l) = 1 

verificando-se que o seu intervalo de definição não é limitado superiormente. Finalmente para 
qualquer k > 0, defina-se 



y(t) se í<l [ ^/12(í- 1) 2 + 1 se í < 1 

n(í)= " 

j/ fc (í) se t > 1 ^ ^(í- 1) 2 + 1 se í>l 

é diferenciável em M, verifica a equação diferencial, verifica a condição inicial e está definida 
em IR. Construimos assim uma infinidade de soluções do PVI definidas em HL 

(iii) Como já referimos, o PVI tem solução única, enquanto (t,y(t)) não atinge a fronteira 
de D, isto é enquanto t < 1. No entanto quando í = l o teorema deixa de ser aplicável pois / 
não verifica as suas hipóteses. Numa vizinhança do instante t = 1 deixamos de poder concluir 
algo a partir do Teorema de Picard, visto que se escrevermos o nosso problema na forma 

y = f(t,y) 

2/(1/2) = 2 

temos que: 

\-y 2 

e a função não está definida em t = 1. 



5. Mostre que o seguinte problema de valor inicial: 



.3/(0) = 1, 

tem uma única solução y(t), definida para t G [0, +oo[, e calcule lim y(t) 

t— >+oo 



Sugestão: Não tente resolver a equação diferencial. Considere a função u{t) definida por 




Uma vez determinada a função u(t), mostre que 



dy > i 

dt ' 



3( U (Í)) 2 +^(Í+1) 2 



e integre esta relação entre Oeí. 

Resolução: 

Definindo 

o domínio de / é 



f(t,y) = 



3 (y(t)f + i/(t + iy 



D={(t,y)eR 2 : 3 (u(t)) 2 + {j {t + l) 2 + o} = R 2 \ {(í, y) = (-1, 0)} 

Comecemos por mostrar existência e unicidade de solução local. Verifica-se facilmente que 
tanto / como df /dy são contínuas em D e visto (to,yo) = (0, 1) G D, o Teorema de Picard 
assegura existência de uma solução do PVI y(t), definida para t numa vizinhança, de to = 0, 
isto é, a solução y(t) existe e é única para t G / =]a,/3[. Visto que G /. podemos concluir 
que y(t) existe e é única para t G / = [0,/?[. Falta mostrar que /3 = oo, isto é, que nem (í, y(í)) 
atingem a fronteira de D para qualquer í > 0, nem \y(t)\ —> oo em tempo finito. Como não 
conhecemos a solução do PVI, teremos que usar um teste de comparação. 
Por um lado 

f(t,y)>0 , V(t,y)ED (2) 

Considere-se o PVI 

( ú = 

\ u(0) = 1 
A sua única solução é u{t) = 1, e como consequência de (2) 

y(t) > 1 , Vi G M (3) 

De (3) podemos concluir que 



y(t) y^ O para todo t > 0, e 

y(t) é limtada inferiormente em [0, oo[. 
Temos ainda que mostrar que y(t) "não explode" em [0, oo[. Visto 

y 2 > , Vy g M 
podemos escrever que 



(í+ l) 2 /3 

Considere-se o PVI 



f(t,y)< T , ,\mk > Ví>0 , yGR (4) 



(t + l) 2 / 3 



v(0) = 1 

A sua única solução é v(t) = 3\/t -\- 1 , e como consequência de (4) 

y(t) < 3\/t + l , Vi G R (5) 

Dado que a função v(t) está definida em [0, oo[ (na verdade em ] — l,oc[ mas para a nossa 
análise só precisamos de saber o que se passa para t > 0), podemos então afirmar que y(t) "não 
explode" no intervalo [0, oo[. Conclui-se que y(t) está definida para t G [0, oo[. Finalmente para 
concluir algo sobre o seu limite quando t —>■ oo, temos 

lim u(t) < lim y(t) < lim v(t) 

t— too t— >oo t— »oo 

ou seja 

1 < lim y(í) < oo 

t— *oo 

Estas desigualdades não nos permitem concluir algo sobre lim í ^ 00 y(í). Para estimarmos este 
limite de modo mais preciso, necessitamos de minorar f(t,y) de forma menos grosseira. Para 
tal, usando (5), podemos concluir que 

/(*. y) > — 1 > — ; 1 (6) 

' - 27{/(7TT)2 + (í + 1)2/ 3 - 28</(í + l) 2 

Considere-se o PVI 

1 

w 



28(í+ l) 2 /3 

w(0) = 1 
A sua única solução é w(t) = ^v^í + 1, e como consequência de (6) 



Z/(í) > — ^í + l , Vi G M (7) 



Finalmente (usando (5) e (7)), concluimos que 

3 



logo: 



y/t + l < y(t) < 3\/t + l , Ví > 

28 



lim y(t) = +oo 

t— S-OO 



